Slow dynamics in a turbulent von Karman swirling flow 



A. de la Torrc^'B and J. Burguete^'0 

' Departamento de Fisica y Matemdtica Aplicada 
Universidad de Navarra, P.O.Box 177, E-31080 Pamplona, Spain 

We present an experimental study of a turbulent von Karman flow produced in a cylindrical 
container using two propellers. The mean flow is stationary up to Re = 10*, where a bifurcation 
takes place. The new regime breaks some symmetries of the problem, and is time-dependent. The 
axisymmetry is broken by the presence of equatorial vortices with a precession movement, being the 
velocity of the vortices proportional to the Reynolds number. The reflection symmetry through the 
equatorial plane is broken, and the shear layer of the mean flow appears displaced from the equator. 
These two facts appear simultaneously. In the exact counterrotating case, a bistable regime appears 
between both mirrored solutions and spontaneous reversals of the azimuthal velocity are registered. 
This evolution can be explained using a three-well potential model with additive noise. A regime of 
forced periodic response is observed when a very weak input signal is applied. 

PACS numbers: 05.45.-a, 47.27.-i, 47.32.-y, 47.65.-d 



Introduction - Turbulent flows arc ubiquitous in na- 
ture: ranging from small scales (heart valves, turbulent 
mixing) to very large scales (clouds, tornadoes, oceans, 
earth mantle, the sun and other astrophysical problems), 
turbulence is present in many applied and fundamental 
physical problems But, in spite of the attention it 
has received, there are still many open questions: the 
emergence of coherent structures, as vortices, in fully de- 
veloped turbulence or the rise of different bifurcations on 
the mean fiowQ. 

Here we will analyze a particular configuration, the von 
Karman swirling flow, where two different propellers are 
rotating inside a cylindrical cavity. These ffows have been 
studied analytically 0, 01 , numerically and exper- 

imentally d, 0, ■ Recently it has been shown that they 
can present multistability and memory effects Thus, 
this configuration is a natural candidate to study the in- 
stabilities that appear in fully developed turbulence and 
the role played by symmetry breaking and coherent struc- 
tures. They have been also used in magnetohydrodynam- 
ics (MHD) experiments looking for the dynamo action 
with recent successful results Il2l and with a very rich dy- 
namics of the magnetic field [l3|. Because of turbulence, 
the whole MHD problem can not be studied numerically: 
a usual approach is to deal with stationary mean flows in 
the so called kinematic dynamo scheme. Although these 
numerical studies have been used to predict thresholds of 
the dynamo action[3,5E,[3l, real flows can present slow 
dynamics compared to the magnetic diffusion time that 
can not be neglected and should be taken into account 
in the numerical codes. 

In this paper we will focus on the slow dynamics that 
appear for very large Re numbers. The flow presents 
a non-trivial alternation between two states that break 
symmetries of the problem. This dynamics, which can 
be assimilated to a Langevin system with a classical ex- 
ponential escape time and forced periodic response. 

Experimental setup - The experimental volume (FiglT|) 
is a closed, horizontal cylinder whose diameter D = 2R — 
20 cm is fixed while the height H can be modified con- 
tinuously. Two propellers are placed at both ends, with 



radius Rprop = 8.75 cm and 10 curved blades, each blade 
with a height of 2 cm and a curvature radius of 5.0 cm. 
Using the standard cylindrical coordinate system, the 
north N (resp. south S) propeller placed at z = H/2 
(resp. z — — iJ/2) has negative (resp. positive) azimuthal 
velocity. The propellers are impelled by two indcpen- 





FIG. 1: Experimental setup, (a) Horizontal cylinder with 
the propellers, inside the tank. The north propeller is at 
z — H/2 and the south is at z = —H/2. (b) Photograph of 
the propeller, (c) Scheme of the south propeller viewed from 
the equatorial plane. The rotation sense with the convex side 
sets the azimuthal velocity as positive. 

dent motors of 3 kW total power, allowing a rotation 
frequency in the range / = — 20 Hz. The frequency of 
the motors is controlled with a waveform function gener- 
ator and a PID servo loop control. The cylinder is placed 
inside a tank of 150 1 of volume in order to avoid optical 
problems and to assure the temperature stability. The 
fluid used is water at 21° C. 

The measurement of the velocity field is performed us- 
ing a LDV system (with a chosen spatial resolution of 
1 cm and temporal resolution up to 100 kHz) and a PIV 
system (spatial resolution of 1 mm and temporal reso- 
lution of 15 Hz). The LDV system allows the measure- 
ment of two components of the velocity field (axial Vz 
and azimuthal vg), while the radial component (vr) is 
obtained by mass conservation. The velocity field ob- 
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FIG. 2: (Color online) (a) Stream-vectors Vr,Vz in the plane 6 = 0. (b) Contourplot of Ve. The north propeller is on the right, 
the south propeller on the left, (c) Vortex visualized with air bubbles, (d) Gray zones indicating the regions where the figures 
(a,b,c) were obtained in the cell. 



tained in this way is consistent with the PIV measure- 
ments in an axial, horizontal plane {9 = 0). Both tech- 
niques are based on the displacement of small particles 
(14//m, p = 1.65 g/cm'^) inside the bulk of the fluid. 

In our experimental setup three parameters can be 
modified. The first one is the Reynolds number Re = 
RVprop/v, defined using the propeller's rim velocity 
Vprop = '^T^Rpropf ■ This number can be varied contin- 
uously in the range Re ~ 10"^ ~ 10^. The second one is 
the aspect ratio of the experimental volume, T = H/D, 
set to 1 in this experiment. The last one is the frequency 
disymmetry, defined as A = {Jn - fs)/ Un + Is), where 
/jv and fs are the frequencies of each propeller, north and 
south. This parameter is varied in the range | A |< 0.1. 

Results - For the Re range explored, the flow is in 
fully-developed turbulence regime. The mean flow V = 
{v) = {Vr,Vg,Vz) represented in Figl2{a,b) is obtained 
with the LDV system, averaging velocity series longer 
than 300 times the period of the propeller. The mea- 
surement is done in the plane 9 = with Re = 3 10^ 
and A = 0. The flow is divided into two toroidal cells, 
each of them following its propeller (positive azimuthal 
velocity in the south, negative in the north). In each 
cell, the flow is aspired trough the axis towards the pro- 
pellers, where it is ejected to the walls. The flow then 
returns along the cylinder's wall, approaching the axis 
near the equatorial plane. Using PIV measurements, the 
instantaneous field can be obtained, and no traces of the 
mean flow are observed. This is due to the high turbu- 
lence rate {rms value over the mean value) which vary 
between 60 — 150%, depending on the spatial position 
and the velocity component measured. 

This mean flow does not preserve the symmetry around 
the equator (a n rotation around any axis in the z = 
plane, i.e. Rj^ symmetry). In the case presented in 
figure [2Ua,b) the cell near the north propeller is bigger 
than the south one. The broken symmetry is recovered 
when the mirror state (south cell bigger than the north 
cell) is considered. Each state (labeled as 'north' N or 



'south' S depending on which cell is the dominant one, 
i.e. bigger) are equally accessible when the system starts 
from rest. This disymmetry is in contradiction with other 
worksQ in which the R-^ symmetry is observed (i.e., the 
frontier between the two rolls is always in the plane z = 
0) , probably due to the presence of baffles or inner rings, 
which would enforce this condition. Other experiments 
at much slower Reynolds numbers [3] (R^ £ 600) have 
shown that the axisymmetry can be broken producing 
near heteroclinic orbits, but preserving the R^^ symmetry. 

When larger tracers (air bubbles) are used, coher- 
ent structures as vortices (Figl^t) are visualized. These 
structures mostly relay in the dominant cell, so the az- 
imuthal velocity of the vortices arc directly related with 
the azimuthal velocity of the dominant cell. These vor- 
tices have a characteristic size D^ortex ~ 5 cm and they 
appear simultaneously to the disymmetry of the mean 
flow. Previous works in similar configurations 0] showed 
the formation of a static vortex in the equatorial plane 
(z = 0) for much lower Re, inaccessible with the present 
configuration. 

The state iV or 5 of the system can be characterized 
using different variables. One is the position of the fron- 
tier zq, defined as the z position in where Ve = (in 
Figl2j3, zq ^ —20 mm). Another possibility is measuring 
Vg'^, the mean azimuthal velocity at an equatorial point 
near the wall (r = 0.9R, z = 0) with the LDV system (in 
FiglSbj Vg"^ ^ —0.4 m/s). This mean velocity is stable in 
time and proportional to the propeller's rim velocity. For 
Re < 10* we find that the normalized azimuthal velocity 
(Ug = VP /Vprop) is nearly null so the mean flow is al- 
most symmetric. As the Re is increased the disymmetry 
becomes more notorious, until a plateau {Ug = const) is 
reached for Re > 10^. 

For this range of Re the system can spontaneously 
jump from one state to the other {inversions) . In Fig[3l^a) 
the instantaneous azimuthal velocity of the equator is 
plotted versus time {tacq ~ 5 10"* s 14 h) for Re = 
2.7 10^, A = 0. The typical transition time is about 10s, 
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FIG. 3: (a) Inversions of the azimuthal velocity at a point 
near the wall in z = for A = 0, 7?e = 2.5 10^. Sohd white 
line: filtered data using a low-pass filter with a cutoff fre- 
quency fcut = 0.025 Ilz.(b) PDF of ue for each state {N or 
S) obtained for A = and various Re. 



while the time between inversions can vary from min- 
utes to hours. These time scales are much slower than 
the period of the propeller {tacq 



2.4 10^ Tprop, with 



Tprop — 1/ fprop — 0-2 s). 

The probability density function (PDF) of these states 
can be computed breaking up the data series in time 
intervals: a low pass filter is applied (figlSja), white 
line) that allows to differentiate between N and S and 
split the signal. The shape of these PDFs does not de- 
pend on the cutoff frequency except for extreme values 
{fcut ^ or fcut — > samphng rate). FigEKb) shows the 
PDF of the normalized instantaneous azimuthal velocity 
{ug = vg/Vprop) for the two states: north (with negative 
most probable ug) and south (with positive most proba- 
ble velocity). Each distribution [pN or ps) is described 
as the superposition of two gaussians: 



PN,siue) = Go + Gn.s = 



A 



N,S 



■ exp 




exp 



2a2 



UN,S) 



(1) 



with Aq + Apf^s = 1- 

Both distributions pn and ps share the same proper- 
ties: one of the gaussians Go has zero mean, while the 
other Gjv,s is centered around a finite value (ImatI = 
I'^sl 7^ 0) that increases slightly with the Re number. 
The amplitudes Aq and Aj^ s have a stronger depen- 
dency on the Re (Fig[3^). For low Re the PDF arc nearly 
gaussians (Aq ~ 1) while for large Re (plateau) the non- 
symmetric gaussian becomes dominant A^^s 3> 7^ 0. 
The zero mean gaussian Gq is due to residues of the sym- 
metric flow and the other gaussian Gn,s is related to the 
displacement of the vortices around the equator. 
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FIG. 4: (Color online) (a) Experimental amplitude of the 
gaussians Ao (x,- -) and An,s (+, — ) in Eq[T] us. Re. The 
lines are only plotted to indicate the trend and are not ob- 
tained from any fitting, (b) PDF of ue for the S state for 
the experimental data series Re = 2.5 lO'' of Fig|3] (solid line) 
and numerical data (dashed line) using the model (eq[5)) with 
e = —0.05, g = 2, B = 1. (c) Numerical amplitudes of the 
uo ( — ) and un,s{- -) fs. e. Each point corresponds to an 
ensemble average of 100 realizations. The noise intensity in- 
creases linearly with e, in such a way that 5 = (resp. 1) 
when e — —1 (resp. —0.05). (d) PDF of the escape times for 
experimental data of Fig. [3] (+) and for the numerical data 
presented in (b) (dashed line, time unit Taim = 391s). 



According to this description, the system visit three 
different regions in phase space around ug = [0,un,us]. 
A simple model based on a three well potential (one for 
the symmetric case = and the two others for the 
asymmetric states u^ g) will describe this dynamics: 



iig = eug + g Ug 



(2) 



where ^(t) is a noise distribution with noise level B (play- 
ing the role of the turbulence rate) and g controls the 
relative depth of the potential wells. Thus, the N state 
(resp. S) will appear when the system is wandering be- 
tween the wells uq and un (resp. uq and ug). The pa- 
rameter e is varied in the range (— g^/4 < e < 0) where 
the three solutions ug = [0, ±{g/2+ {g'^/A + e)0-5)0-5] are 
stable. 

Different runs using an Euler-Maruyama scheme [l3| 
were performed in order to recover the dynamics: for 
small e the dynamics is confined to the region around 
Ug = Q whereas for e ^ the numerical evolution 
presents spontaneous inversions. In this later case the 
PDF of each state can be computed and compared to ex- 
periments (fig. Hlb). The characteristic doubly bumped 
distribution is obtained, but in the numerical distribu- 
tion the queues arc not symmetric due to the shape of 
the Q-th order potential in the neighborhood of us- The 
relative weight of each one of the solutions in the numer- 
ical PDFs can be calculated (fig. [He) and compared with 
the experimental amplitudes of fig. |3)a. 
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The distribution of the times the system stays in one 
state {residence times) follows an exponential decay law 
(Kramer's escape rate [l8, 1^, 2^): 



p(t) = l/ToexpH/ro) 



(3) 



where Tq is related to the intensity of noise. In figure IHd 
we present the experimental residence time for the data of 
fig.[3]a. The experimental data have a characteristic time 
To = 1484s = 7020 Tp^p = 0.19 r^, being = R^/iy the 
diffusion time scale. 

As a consequence of the dynamics, a natural ques- 
tion that arises in this problem is the response of 
the experiment to an external forcing. When a sinu- 
soidal modulation is applied to the frequency of one of 
the propellers, the rim velocity evolves as Vprop{t) = 
Vpj.gp {1 + 2Ao cos{ujt)) where Aq is the maximum fre- 
quency disymmctry. 
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FIG. 5: (a) Position of the frontier {zq) varying the difference 
between the propellers' frequency (A). The system presents 
a bistability at A = 0. (b) Temporal evolution of the instan- 
taneous velocity of an equatorial point near the wall, with 
A = 0.026 cos(0.027rt). Re = 2.610^. The system remains 
slaved to the input signal. The time scale is two orders of 
magnitude smaller than that of figure [S] a 

A static forcing (w = 0, Aq ^ 0) was applied to obtain 
the minimum amplitude needed to induce an inversion 
in the experiment. In Fig[5Ua) zq is represented vs A. 
A disimmetry of only A ~ 0.0025 is needed to make the 
system jump between both states. No inversions have 
been observed with A 7^ 0, so up to our precision and 
measurement times no hysteresis is found, and the system 
only shows a bistability in A = 0. Further investigation 
will determine if Tq < cx) when A 7^ 0. 

When the harmonic forcing is induced, preliminary re- 
sults show how for low / only a small amplitude in the 
input signal is needed to obtain the resonance of the sys- 
tem; Periodic inversions are observed with the same fre- 
quency of the forcing. Figure [3lb shows how the system 
is slaved to a signal with Ao = 0.026 and / = 10~^ Hz. 
In the range 10~^ > ///prop > 10^'^ it is necessary to 



increase A to observe the synchronism. For higher /, 
the system can not follow the forcing: the inertia of the 
flow limits the response time of the inversion Tres '^-^ 20 s 
~ 2.6 10"'^ r^. Further work experimental and mmrerical 
is under run in this direction. 

Conclusions - We have presented experimental evi- 
dence of a bifurcation in a turbulent system that breaks 
symmetries and produces slow dynamics. Two possible 
mirror states are equally accessible, with random (natu- 
ral) or periodic (induced) inversions. The slow dynamics 
can be characterized using a very simple Langevin model 
in spite of being in a fully developed turbulence regime. 
Finally, when a very weak amplitude input modulation is 
applied to the system, a forced periodic response appears 
in the turbulent flow. One open question is either these 
natural inversions could be present in dynamo experi- 
ments for very long temporal series (Tg ~ 1000s) when 
A = 0. 
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